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A typical signature of charge extraction in disordered organic systems is dispersive transport, which implies 
a distribution of charge carrier mobilities that negatively impact on device performance. Dispersive 
transport has been commonly understood to originate from a time- dependent mobility of hot charge 
carriers that reduces as excess energy is lost during relaxation in the density of states. In contrast, we show via 
photon energy, electric field and film thickness independence of carrier mobilities that the dispersive 
photocurrent in organic solar cells originates not from the loss of excess energy during hot carrier 
thermalization, but rather from the loss of carrier density to trap states during transport. Our results 
emphasize that further efforts should be directed to minimizing the density of trap states, rather than 
controlling energetic relaxation of hot carriers within the density of states. 



While natural photosynthesis transfers electrons through a cascade of energy states, artificial photovol- 
taic systems must extract photogenerated charges to the electrodes, and despite recent performance 
gains 1 , fundamental questions about this charge extraction still remain unanswered. There has been 
intense scrutiny of the mechanisms of charge generation and the impact of above-bandgap photon energy 2 " 4 , 
however, this level of attention has not extended to studies of the extraction of such 'hot' charge carriers, despite 
the fact that efficient charge extraction is crucial for device performance 5 . 

The most characteristic feature of charge transport in disordered systems is the dispersion of the charge carrier 
movement velocities 6 . Dispersive transport harms device performance because the slowest carriers bring down 
the average mobility 7 , and consequently, the vast majority of novel organic semiconductors remain inapplicable 
for efficient devices. Moreover, the detrimental effects of dispersion are exacerbated by the inhomogeneities in 
film thicknesses caused by the targeted low cost deposition methodologies, because the transit time distributions 
become dramatically longer and more dispersed in regions of increased thickness. 

Dispersive transport in organic semiconductors is usually thought to be caused by the energetic relaxation of 
hot charge carriers within their density of states 8 . Spectroscopic measurements and Monte Carlo simulations have 
revealed energetic relaxation extending even to the microsecond timescales, where it could be relevant to bulk 
charge transport 910 . Even if the bulk of the energetic relaxation were to occur on very fast timescales, there is still 
the question of whether residual thermalization might continue to long, microsecond timescales. This energetic 
relaxation is often understood to cause a time- dependent mobility and therefore explain dispersive current 
transients 1112 , yet we will show here that this commonly-used model is inconsistent with our observations in 
high efficiency organic solar cell materials. Instead, there is an alternative mechanism for the creation of a 
distribution of carrier velocities, namely, via trapping. This observation has a very direct impact on the numerous 
models, theories and experimental results describing dispersive charge transport in disordered organic semi- 
conductors. Furthermore, it points to a new strategy for improving charge transport "management" in devices 
such as organic solar cells. 

The classic signature of dispersive transport is a time-of-flight photocurrent signal that decays with time even 
before the carriers have transited through the film 13 . This decay in photocurrent can occur due to two mechan- 
isms, a reduction in carrier mobility, and/or a reduction in the number (or concentration) of moving carriers. The 
former, a time- dependent hot carrier mobility, is presently commonly believed to be the cause of dispersion in 
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organic semiconductors 81112 and it is usually understood to originate 
from a loss of energy as carriers thermalize within their density of 
states 14 " 16 . Higher energy carriers are expected to have a higher hop- 
ping probability, and hence a higher velocity 17 " 19 , so the thermaliza- 
tion within the density of states causes the carrier mobility to decline. 
Recent studies 910 have reported mobility thermalization times on the 
order of microseconds. However, an alternative explanation for the 
decaying transient photocurrent, which is less commonly accepted in 
organic semiconductors, is a time dependent concentration that can 
arise if carriers are gradually lost to traps 20 " 22 . The photocurrent 
signal will continue to reduce as long as the net concentration of 
moving carriers continues to decrease. If that physical process pre- 
vails, there can be decaying photocurrent despite the moving carriers 
having a constant drift velocity. Additionally, if the cause of disper- 
sion is trapping, then it will influence all devices, even those which 
operate in the dark 23,24 . 

In this article, we demonstrate that a time-dependent hot carrier 
mobility cannot explain the dispersive transport in our studied bulk 
hetero junction solar cells. We address this issue by performing tran- 
sient photoconductivity experiments in which we vary the transit 
time by changing the electric field and/or device thickness. The 
expectation is that if the dominant cause of dispersive transport is 
mobility relaxation, then the average mobility and the amount of 
dispersion should vary with the electric field and/or film thickness, 
because longer transit times will allow for more relaxation to occur. 
Conversely, if the dominant effect is trapping, then it is the concen- 
tration of carriers which is changing in time rather than their mobil- 
ity, and consequently, the average mobility and the dispersion range 
should not vary with film thickness or electric field. This transit time 
dependence allows these two dispersive mechanisms to be experi- 
mentally distinguished. 

Results 

Numerical Simulations of Resistance dependent Photo Voltage 
(RPV) measurements. Our experiments were made possible by 
the development of a new transient photoconductivity technique 
that we call Resistance dependent Photo Voltage (RPV), which is 
described here and in the Methods section. The experimental 
measurement circuit for RPV is shown in Figure 1. This setup is 
similar to time-of-flight, where charge carriers are photogenerated 
by a short low- intensity laser pulse. A low light intensity is necessary 
so that the electric field inside the device is undisturbed. The 
transient photosignal is determined by the competition between 
two simultaneous processes: the transport of charge carriers inside 
the film, and the response of the external RC circuit. Unique to the 
RPV approach, and in contrast with time-of-flight, the entire 
measurement is repeated at many different load resistances 
spanning the range from differential mode (small R) to integral 
mode (large R). The resistance is varied for two reasons: firstly, to 
visualize the transit times, as will be shown below; and secondly, to 
reveal the slower carrier mobility by amplifying the slower carrier's 
conduction current. The slower carriers produce a much smaller 
current than the faster carriers, and their transit would be buried 
in the noise at resistances that are optimized for the faster carriers. 
Slower carriers have much longer transit times, allowing the use of 
larger resistances, and consequently allowing for their weaker 
electrical signal to be amplified. In this way, RPV bridges the gap 
between differential mode and integral mode time-of-flight, and 
allows measurement of the transport of both types of charge carriers. 

The combination of the RC circuit dynamics, dispersive transport, 
and optical interference effects prevent analytic analysis of the tran- 
sients. To study highly dispersive systems, such as organic solar cells, 
the simultaneous impact of all these effects must be understood. We 
applied numerical simulations to develop this understanding. The 
simulations are described in the Methods section and in the 
Supplementary Information. Typical simulated transients for an 
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Figure 1 | Resistance dependent Photo Voltage (RPV) measurement 
circuit (top) and timing diagram (bottom). A low light intensity 
nanosecond laser pulse is used to photogenerate charge carriers inside (for 
example) the semiconductor junction of an organic solar cell. Low light 
intensity is critical in the RPV experiment to ensure operation within the 
"small charge extraction mode" where the internal electric field 
distribution in the film is not altered by transported charges. After 
photogeneration, the charge carrier transport through the film is driven by 
the built-in or the applied external electric field, and the resulting transient 
photosignal is recorded by an oscilloscope. The transient photosignals are 
measured at various load resistances i?i oa d- 

organic solar cell with dispersive transport are shown in Figure 2 
(a). The transients show two distinct extraction "shoulders," as indi- 
cated by the arrows. The transients at different resistances assist in 
visually identifying the location of these "shoulders." The positions 
of these arrows correspond to the mean transit times required for the 
faster and slower carriers to cross the entire thickness of the film. In 
this simulation, carriers are repeatedly trapped and de-trapped, cre- 
ating dispersion because the total time spent in traps is different for 
different carriers. The resulting distribution of transit times is shown 
at the top of Figure 2 (a), and its approximate width is indicated by 
the shaded background. It can be seen that the RPV technique allows 
the mean charge carrier mobility to be obtained even in the presence 
of strong dispersion. 

In addition to shallow traps that cause dispersion, we also con- 
sidered deep traps that immobilize carriers for times much longer 
than the transit time of either carrier. Long lived trapping is typical in 
disordered organic semiconductors 22,25 , because many organic mate- 
rials behave as unipolar conductors, and solar cells often have 
strongly imbalanced mobilities 26 . In these cases, repeated photoge- 
neration adds more trapped charge in the form of the immobilized 
charge carriers, which might accumulate with every repetitive laser 
shot, redistributing the electric field and distorting the measurement. 
Figure 2 (b) shows simulations of this film charging for the case of 
fast Langevin-type recombination under repeated laser shots, as 
would arise from the presence of deep trap states far inside the 
forbidden energy gap 27 . These are large resistance transients, in other 
words, the measurement circuit has integrated the photocurrent such 
that the peak voltage is proportional to the extracted charge. If the 
extracted charge is decreasing and the extraction time remains con- 
stant, then carriers must be lost to recombination and not due to field 
screening, and hence we conclude that the trapped charges act as 
recombination sites for the mobile carriers. However, the mobility of 
the charge carriers can be determined independently of the trapping 
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Figure 2 | Numerically simulated RPV transients in the case of (a) dispersive transport caused by shallow traps, and (b) film charging caused by deep 
traps, (a) In the case of dispersive transport, the extraction "shoulders" approximately correspond to the mean charge carrier mobility, (b) In the case of 
deep traps, the film becomes charged and the magnitude of the RPV transient is reduced in subsequent shots of the laser, but the transit "shoulder" 
remains unhindered which allows for reliable charge carrier mobility estimation. 



effects, because the rapid Langevin recombination prevents the 
build-up of large amounts of trapped charge that would disturb the 
transit time. 

Experimental measurements. We chose to study the well-known 
photovoltaic blend 28 of poly[AT-9"-hepta-decanyl-2,7-carbazole-a^- 
5,5-(4',7'-di-2-thienyl-2M',3'-benzothiadiazole)] (PCDTBT) and 
[6,6]-phenyl-C 70 -butyric acid methyl ester (PC 71 BM) in an 
optimized blend ratio of 1:4 by weight 29 . This blend is ideally 
suited to this study because its amorphous nature allows the 
elimination of any film thickness dependent morphology 30 . In 
order to see the generality of observed effects, we have also done 
the same experiments on poly[[4,8-bis[(2-ethylhexyl)oxy]benzo 
[l,2-b:4,5-b']dithiophene-2,6-diyl] [3-fluoro-2-[(2-ethylhexyl)carbonyl] 
thieno[3,4-b]thiophenediyl]] (PTB7):PC 71 BM blends, the results of 
which are shown in the Supplementary Information. The thin film 
(active layer thickness of 75 nm) PCDTBT:PC 71 BM solar cell 
exhibited a power conversion efficiency of 6.3% under standard 
AM1.5G illumination, while the PTB7:PC 71 BM blends reached 
7.7%. Current-voltage curves for both devices are shown in 



Supplementary Figure 5. None of the optimized PCDTBT or PTB7 
based devices demonstrated any significant film morphology 
inconsistencies in the range of studied film thickness. (See the 
Methods section for the details of the fabrication and the 
Supplementary Information for characterization of photovoltaic 
performance). The presence of dispersive transport was confirmed 
by time-of-flight experiments on thick films (Supplementary Figure 
6). No photocurrent plateaus were observed; the transients decrease 
with time as is typical of dispersive systems. 

Figure 3 shows the recorded RPV transient signals for a 
PCDTBT:PC 71 BM solar cell. All transients were recorded at near 
to short-circuit conditions. This remains true even at large resis- 
tances, because the maximum photovoltage occurring during the 
transient is substantially less than the built-in voltage. The first 
shoulder marks the arrival time of faster carriers (27 ns), which is 
attributed to electron transport since the time scale is similar to that 
measured for PC 7 iBM (please refer to the Supplementary 
Information for measurements on PC 71 BM). The second shoulder 
is less well defined due to the strongly dispersive nature of hole 
transport in this system, but marks the arrival of the slower carriers 
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Figure 3 | Experimentally measured RPV transient photo-signals in an 
optimized PCDTBT:PC 71 BM solar cell. Mean electron (faster) and hole 
(slower) transit times are marked, from which the respective mean 
mobilities are estimated. The dispersive nature of charge transport in the 
studied solar cells is highlighted by shaded boxes marking the range of 
carrier arrival times. Thin curves show recorded data, while bold lines show 
data smoothed by adjacent averaging. The short timescales for large 
resistances were omitted for clarity. 

(2.59 |is). Mean electron and hole mobilities were determined from 
the shoulders in the transients, as indicated by arrows in Figure 3, 
with the approximate spread of arrival times indicated by the shaded 
boxes (corresponding to the regions where the transients deviate 
from the dotted lines). The edges of these shaded boxes give the 
"fastest" and "slowest" case transit times, from which we obtained 
the dispersion range in the mobilities for each species. This range is 
an essential feature of the dispersive transport exhibited by this sys- 
tem, because a single mobility value does not correctly quantify the 
transport when the system is dispersive. We measured the mean 
electron mobility to be 2.9 X 10~ 3 cm 2 V" 1 s _1 with a dispersion 
range from 1.1 X 10" 3 cm 2 V" 1 s" 1 to 4.5 X 10" 3 cm 2 V" 1 s" 1 and the 
mean hole mobility to be 3 X 10" 5 cm 2 V" 1 s" 1 with a dispersion 
range from 9.2 X 10" 6 cm 2 V" 1 s" 1 to 7 AX 10" 5 cm 2 V" 1 s" 1 . Despite 
the high level of dispersion observed here (the hole dispersion range 
covers nearly an order of magnitude), the OPV device still maintains 
good performance. However, further work is necessary to identify the 
impact of the dispersion range on the performance of solar cells. 

Next, we studied the impact of photon energy on the hot charge 
carrier transport, because any relaxation effects are likely to be 
dependent upon the initial energy. This is important because of 
recent suggestions that excess above-bandgap energy may assist exci- 
tonic dissociation 4 , although the methodology of that observation 
has been challenged 31 . We note that quantum yields have been shown 
to be independent of the energy level of the excited state, suggesting 
that hot excitons are indeed not beneficial for exciton separation 32 . 
Nevertheless, hot charge carriers - rather than excitons - might also 
possess excess energy and shape the internal quantum efficiency 
spectra; therefore, it is important to clarify these effects, aiming for 
improvement in the charge extraction of typical low mobility organic 
materials. In the past the absence of hot charge carrier effects has 
been observed indirectly 33 . Numerical simulations predict that RPV 
is independent of optical interference effects (Supplementary Figures 
2 and 3), allowing direct and unambiguous measurement of any hot 
charge carrier effects that may be present. RPV transients were mea- 
sured at two different photon energies, 3.49 eV (355 nm) and 
2.33 eV(532 nm). The results are plotted in Figure 4, showing nearly 
identical transients resulting from laser excitation at the two different 
wavelengths. The photon energy independent mobility suggests that 
excess energy plays a minimal role in dispersive transport, since 
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Figure 4 | RPV transients measured on a 75 nm PCDTBT:PC 71 BM solar 
cell using two different laser wavelengths: 355 nm (3.49 eV) and 532 nm 
(2.33 eV). The nearly identical transient responses directly demonstrate 
the absence of hot carrier effects in this system. 



carrier thermalization (if it is present) must happen in time scales 
much shorter than the transit time. 

To further confirm that the dispersion in hot carrier mobilities is 
not caused by the thermalization of carriers, we studied the electric 
field and film thickness dependence. Longer transit times should 
allow more time for thermalization, thus influencing the result if 
the dispersion is due to carrier relaxation. The results are shown in 
Figure 5; the Supplementary Information includes a selection of the 
transients from which these mobilities were estimated. The mobili- 
ties and dispersion ranges are completely independent of electric 
field and photon energy [Figure 5 (a)], suggesting that trapping 
mechanisms are more significant than relaxation mechanisms. The 
lack of electric field dependence is in contrast with the Poole-Frenkel 
dependence reported in pristine PCDTBT 34 . This is an unexpected 
result, because in disordered organic systems significant electric field 
dependence is typically observed, even at relatively low values of 
electric fields 34 , which is thought to originate from hopping-type 
charge transport. Further studies of the temperature dependence, 
and measurements on other systems, have to be performed in order 
to clarify the origin of this observation. Additionally, we observe that 
the mean mobilities and dispersion ranges are nearly independent of 
the film thickness [Figure 5 (b)]. We attribute the small changes in 
mobility to device-to-device variations that result from the fabrica- 
tion process. The thickness independence of the mean mobilities and 
dispersion ranges further support the claim that the dispersion is 
caused by traps instead of relaxation. A charge carrier density 
dependence in the mobility even at low concentrations has been 
observed in P3HT:PCBM blends 24 , and we note that a concentration 
dependence might cause dispersion as carriers gradually become 
trapped and the density decreases. We do not exclude the possibility 
of a density dependence here. However, in our measurements, 
increasing thickness corresponds to lower densities because the 
amount of photogenerated charge was always less than CU, which 
is inversely proportional to thickness. Consequently, the thickness 
independence in the mobility implies that there is negligible density 
dependence at the concentrations probed here. 

Further measurements were also performed on solar cells made 
with PTB7 blends. The results show the same conclusions as the 
PCDTBT blends: the mean mobility and dispersion ranges are inde- 
pendent of film thickness, applied electric field, and photon energy 
(Supplementary Figures 8, 10, and 11). The results reported here 
appear to be generally applicable and are certainly not specific to 
PCDTBT blends. 
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Figure 5 | Electron and hole mobilities measured in PCDTBT:PC 71 BM 
solar cells. The error bars show the dispersion ranges. Carrier mobilities 
and dispersion ranges are independent of electric field and photon energy 
[panel (a)], and nearly independent of film thickness [panel (b)], 
demonstrating that carrier thermalization cannot account for the 
dispersive transport in this system. Consequently, dispersion is caused by 
trapping. 



Discussion 

Charge transport in the studied operational OPV blends is strongly 
dispersive, as demonstrated by the decaying time-of-flight photo - 
current transients in thick devices (Supplementary Figure 6). These 
time-of-flight transients were recorded in a regime where drift 
dominates over diffusion, so the current density is described by 
j = eE (n/Li n + pfi n ), where e is the charge of an electron, \i n and \i p 
are the electron and hole mobilities, n and p are the carrier concen- 
trations, and E is the electric field. The observation of a decaying 
photocurrent density j can be explained by two mechanisms: ther- 
malization (a time dependent mobility, fi), and/or trapping (a time 
dependent concentration of moving charge carriers, n). These 
mechanisms are schematically illustrated in Figure 6, from which it 
can be seen that either model would result in dispersive photocurrent 
transients. We found no evidence of thermalization-type effects on 
the timescales comparable with those involved in charge transport. 
Figure 4 directly demonstrates that that excess energy of hot carriers 
has essentially no contribution to mobility or dispersion. In Figure 5, 
we demonstrate that the dispersion range is independent of the 
applied electric field and changes very little with thickness. If ther- 
malization on transport time scales 16 were the cause of the dispersion, 
then modifications to the transit time should change the mean 



mobility and/or dispersion range by varying the time available for 
relaxation. Such a variation was not observed, and hence we exclude 
thermalization as the mechanism of the dispersive transport. Any 
relaxation processes must be much faster than charge transport, so 
that the distance covered by charges as they relax is insignificant 
compared with the film thickness, and hence the relaxation has neg- 
ligible contribution to the overall dispersion. With relaxation 
excluded, the only remaining mechanism is a reduction in the con- 
centration of moving carriers, therefore, we conclude that trapping is 
the primary cause of the dispersion in these systems. This challenges 
the widely- used model of hot carrier relaxation within the density of 
states. Consequently, dispersive transport potentially impacts on the 
many different devices that employ films made from disordered 
semiconductors, including those that operate in the dark or at 
steady- state conditions. 

In conclusion, electron and hole mobilities and their dispersion 
ranges were measured simultaneously using the RPV technique in a 
high efficiency narrow optical gap polymer/fullerene system 
(PCDTBT:PC 71 BM). We found that the transport of electrons and 
holes are both strongly dispersive in these thin, efficient solar cells. 
We introduced the dispersion range as a parameter to quantify 
charge transport, since a single mobility value is insufficient to prop- 
erly characterize a dispersive material. We directly observed the 
absence of "hot carrier" effects on time scales relevant to charge 
extraction, and furthermore found that the dispersion is caused by 
trapping rather than thermal relaxation. We have found that the 
widely-used model of hot carrier relaxation within a density of states 
is not the dominant process causing the dispersion in the studied 
solar cells. Furthermore, in contrast with the Poole-Frenkel depend- 
ence previously reported in pristine PCDTBT and other disordered 
systems, the studied solar cell blends exhibit an unexpected negligible 
electric field dependence. While further work is needed to clarify this 
observation, electric field independence may assist in maintaining a 
good fill factor by keeping the mobility higher near the maximum 
power point. The absence of hot carrier effects and an electric field 
independent mobility were also observed in PTB7:PC 71 BM solar 
cells, suggesting that these conclusions may be more generally applic- 
able. This work signifies the importance of localized trap states as 
opposed to thermalization and hot carrier effects in efficient poly- 
mer-based solar cells. Since dispersion arises from trapping, it is also 
important for other types of devices, such as organic field effect 
transistors and diodes. Trap states are relevant whether the carriers 
were injected or photogenerated, and whether the device is in tran- 
sient or equilibrium conditions. Our results suggest that further sci- 
entific research should be directed towards reducing the density of 
trap states rather that utilizing above-bandgap energy for improving 
electronic device performance. 

Methods 

Numerical simulations. The simulations are based on a standard one-dimensional 
drift-diffusion-recombination solver 3536 assuming a negligible amount of equilibrium 
carriers, which is typically the case in organic semiconductors 37 as well as in the 
studied devices. For simulations of dispersive transport, we implemented a multiple 
trapping and release model 20 ' 2138 with an exponential density of localized states. The 
full list of equations are given in the Supplementary Information. 

Solar cell fabrication. 15 Q/sq. Indium tin oxide (80 nm thick, purchased from 
Kintec) coated glass substrates were cleaned in a 100°C water bath with alconox 
(detergent), followed by sonicating in sequence with de-ionized water, acetone and 2- 
propanol for 6 minutes each. Next, a 30 nm layer of poly(3,4- 
ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) was spin-coated at 
5000 rpm for 60 sec onto the cleaned substrates, which were then annealed at 170°C 
for a few minutes in air. For PCDTBT devices, a solution of PCDTBT (purchased 
from SJPC Group) and PC 71 BM (purchased from Nano-C) was prepared by using a 
1 : 4 blend ratio by weight and a total concentration of 25 mg/cm 3 in 
1,2-dichlorobenzene (DCB). Solar cells with four active layer thicknesses, 75 nm, 
230 nm, 270 nm and 410 nm (measured by a DekTek profilometer), were 
fabricated by spin coating. For PTB7 devices, the active layer of PTB7 (1 -Material, 
Mw = 97.5 kDa, PDI = 2.1) and PC 7 iBM (ADS) was prepared as previously 
described 1 resulting in 100 nm, 150 nm, 230 nm, and 700 nm thick films. To 



SCIENTIFIC REPORTS | 4:5695 | DOI: 1 0.1 038/srep05695 



5 



Energy J ^^q 0 




Mobility decreases as carriers relax 
within the density of states. The length 
of the arrow indicates the mobility. 



(a) Thermalisation 



Position 



Energy 



Trap 
states 



Mobility is constant, but carriers are 
progressively lost to trap states 

o \p • ■ 




o o 



(b) Trapping 



Position 



Figure 6 | Schematic illustration of the two pathways to dispersive transport, (a) Thermalisation causes the mobility to decrease with time, whereas (b) 
trapping causes the loss of carrier density. We have shown here that the latter case (trapping) is the dominant effect in the studied solar cells. 



complete the solar cells 1 .2 nm of samarium and 75 nm of aluminium were deposited 
under a 10" 6 mbar vacuum by thermal evaporation. The device areas were 0.2 cm 2 for 
current density versus voltage (J-V) measurements and 3.5 mm 2 for charge transport 
measurements. J-V characteristics were obtained in a 4- wire source sense 
configuration and an illumination mask was used to prevent photocurrent collection 
from outside of the active area. An Abet solar simulator was used as the illumination 
source and provided —100 mW/cm 2 AM1.5G light. 

RPV measurements. A delay/trigger generator (Stanford Research Systems DG535) 
was used to trigger the laser and function generator (Agilent 33250 A) pulses for 
timing control. A pulsed Nd:Yag laser (Brio Quantel) with a pulse length of 5 ns was 
used to generate the carriers. Optical filters were used to reduce the laser intensity for 
the RPV measurements. A function generator was used to apply external voltage 
pulses for electric field dependent mobility measurements. RPV photovoltage signals 
were recorded with an oscilloscope (WaveRunner 6200 A) at various external load 
resistances. RPV transients were smoothed with an adjacent averaging function to 
neutralize the electromagnetic wave oscillations in the measurement circuit. In 
agreement with previous studies done by Clarke et al. 39 , dark-CELIV transient 
responses showed no equilibrium carrier extraction, justifying the application of RPV 
to the studied devices. Optical interference simulations were performed using the 
transfer matrix approach 40 with typical optical constants of PCDTBT/PCBM 
blends 41 . 
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